ON CERTAIN STATISTICAL PROPERTIES OF CONTINUED 
FRACTIONS WITH EVEN AND WITH ODD PARTIAL QUOTIENTS 



FLORIN P. BOCA AND JOSEPH VANDEHEY 

Abstract. We prove results concerning the joint limiting distribution of the renewal time 
of denominators and consecutive digits of random irrational numbers in the case of con- 
tinued fractions with even partial quotients, with odd partial quotients, and for Nakada's 
a-expansions. 



1. Introduction 

Let (a n ) (respectively, (q n )) denote the sequence of digits (resp., denominators of the 
convergents) in the regular continued fraction (RCF) expansion of an irrational number. For 
each R > 1, consider the renewal time ur := min{n : q n > R}, so that q nR -l ^ R < q nR - 
As a consequence of their renewal-type theorem for the natural extension of the Gauss map 
associated with regular continued fractions (RCF), Sinai and Ulcigrai proved the existence of 
the joint limiting distribution of ( gn ^ -1 , a nR ^K, • • • , a-n R +K) with K a fixed nonnegative 
integer [15], as R — > oo. The classical Gauss-Kuzmin statistics give the probability of a 
random x in [0, 1] having a prescribed string of digits in its continued fraction expansion 
at the nth position, for large n; the joint limiting distribution studied in [15, 16] gives the 
probability of a random x in [0, 1] having a prescribed string of digits in its continued fraction 
expansion at the first place where the denominator of the convergent is larger than R, for large 
R. The joint limiting distribution may therefore be considered an analogue of Gauss-Kuzmin 
statistics. Employing an abstract characterization of denominators of successive convergents 
in the regular continued fraction expansion RCF(x) of x, Ustinov succeeded in explicitly 
computing this limiting distribution in the RCF case [16]. 

Sinai and Ulcigrai's result has been subsequently extended to the situation of continued 
fractions with even partial quotients (ECF) by Cellarosi [3]. The ECF limiting distribution 
was further used in the renormalization of theta sums — that is, replacing the theta sum 
2~2 e n ™ with a theta sum of the type 22 e~ wm ' u modulo a rescaling, rotation, and small 
error term — as the map u — > — ^ modulo 2 is closely related to the forward shift of even 
continued fractions. This has lead to some new results about the distribution of normalized 
theta sums and geometrical properties of their associated curlicues [4, 14]. 

This paper studies this type of limiting distributions in the case of three types of continued 
fractions: ECF, OCF (continued fractions with odd partial quotients), and NCF a (the 
Nakada a-expansions, which include NICF, or continued fraction to the nearest integer, as 
a special case). In the ECF case we provide a direct proof of the main result in [3] while 
making the limiting distribution explicit. The analogous problem is also solved in the OCF 
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case, for which no ergodic theoretical approach is known at this time. As in [16], the key 
tool is providing an abstract characterization for pairs of successive convergents in ECF(x) 
and OCF(x), which may be of independent interest. The OCF case is the most intricate, 
because the sequence of denominators of successive convergents in OCF(x) is not necessarily 
increasing as in the RCF, ECF, or NCF a cases. Finally we provide an explicit relation 
between the NCF a limiting joint distribution and the distribution computed in [16]. 

Concretely, for a given type of continued fraction expansion (ECF, OCF, or NCF a ), 
consider the renewal time 

riR = min{n G N : q n > R} = min{n G N : q n ~i ^ R < q n }, R > I, 

and the joint limiting distribution of ( 9n ^~ 1 , w nfl _x, • • • , oj nR+ K) with uJk = (au, e^), for 
fixed K, as R — > oo. Here again, uj^ denote the continued fraction digits and q n denote the 
denominators of the convergents for a given type of CF expansion (see Section 2 for more 
details). 

We will evaluate the Lebesgue measure Cxifxi% 3 ,x4 (R) of the set of numbers x G CI := [0, 1] \ 
Q for which there exist |j,|j7 successive convergents in ECF(x) (respectively in OCF(x)) 
such that for given x\, X2, £3, X4 the following conditions are satisfied: 

R ' ~Q' 2 ' ~Q' X3: ^ ' 

Q'x-P' , Q'x-P' 

^ — - — — ^ x 4 respectively - x 4 ^ — — < 0. (1.2) 
— Qx + P —Qx + P 

In both ECF and OCF situations, we take x±,X2,xs,X4 G (0, l]. 1 In the OCF case, the 
ratio of successive denominators Q/Q' can in fact be any rational number in the interval 
(0, G), but since in the definition of ur we are interested only in Q ^ R < Q' , we can 
restrict our attention to X3 ^ 1 in the definition of C 0,± . The golden ratios G = 1+ 2 V ^ and 
g = ^ = ~ 1 + v/ ^ will be used often. 

The terms gn ^ -1 and in the joint limiting distribution clearly relate to the parameters 
x\ and X2 in the function C. Likewise, the digits UJk in the joint limiting distribution relate 
to the parameters 23 and X4 in C due to equalities (2.4) and (2.6). 

The main result of this paper shows that C E / 0,± (R) has an explicitly computable limiting 
distribution as R — >■ 00. 

Theorem 1.1. The joint distributions C^H%^x^,xJyR) exist as R — > 00 and 

^xi%2,x 3 ,xi{ R ) = + £ (RT 1+£ ), (1.3) 

£°xl% m (.R) = + £ (R- 1 ^), (1.4) 

£xl,X2,X 3 ,X4,(R) = \-0 £ (R 1 l 2+e ), 



1 If any of the parameters equals 0, then C equals as well, so we ignore this degenerate case. 
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where F± = F±(xi, x 2 , X3, X4) and D- L = Di(x\, x 2 , X3, X4) are given by 2 

)U 2 (TXlX 2 X4) ifx 3 ^xix 2 , 
F± = T < (1-5) 
[Li 2 (=Fx 3 x 4 ) - log(l ± x 3 x 4 ) log ^ if x 3 < xix 2 , 

D 2 = F_(xi,x 2 ,x 3 ,x 4 ) - F-(xi,x 2 ,mm{x 3 ,g 2 },x A ), 

£>i = £i+ D Z = Y,IT, with ( L6 ) 

e^i 1^2 

if = f dx [ — — - — - — -, where (1-7) 
Ji/x 2 Jx/(2e+g) y{y ± ^4X) 

X x — 1 



A e = (2£ + g)x 1 , B e (x) = B e ^ X3 (x) = min jx 3 x,xi, 



2f 21 -\ 



The integrals if can be written explicitly as a combination of logarithms and dilogarithms. 

Kraaikamp's metric theory for 5-expansions [6] provides immediate characterizations of 
pairs of successive convergents for such continued fractions, which are obtained from RCF 
only by singularization (see the remark at the end of Section 3 for definition of singulariza- 
tion). In the last section we show how to compute the joint limiting distribution associated 
as above with Nakada's a-expansions [10] for \ ^ a ^ 1. The cases a = 1 and a = \ 
are best known, corresponding to the RCF and NICF (continued fraction to the nearest 
integer). The latter was introduced by Minnigerode [9] and was also studied in [1, 13, 18]. 
Our calculations show explicit connections with Ustinov's RCF distribution. 

2. Basic ECF and OCF properties 
For each x G fl, the ECF (respectively, OCF) expansion of x is given by 

1 

x = = [[(ai, ei),(a 2 ,e 2 ), (03, e 3 ), ■■•]], (2.1) 

ei 

ai H 

e 2 

a 2 H 



e3 

03 H 



where e n G {±1} and all a ra 's are even positive integers (respectively, all a n 's are odd positive 
integers with a n + e n ^ 2). For more details see [5, 6, 8, 11, 12, 13]. As in [5, 8], consider 
the "flipped" continued fraction map Tr> : [0, 1] — > [0, 1] for a subset D of [0, 1], defined by 
T D (0) = 0, T D (l) = 1, and 



T D (x) 



{l/x} if x € (0, 1) \ D, 

l-{l/x} ifxeA 



with auxiliary functions 



jl ifx€[0,l]\A J [l/x] if x G [0, 1] \ D, 

ezHx) = < a D (x) = < 

[-1 ifxGD, ll + [l/x] ]fx€D. 



2 In this paper the convention is that — when a ^ b. 
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Note that 



Consider the sets 



T D {x) = e D (x) 



x 



a D (x)\, Vx 6(0,1). 



°o:= U 

n£2N 



1 



1 



n + 1 n 



D E ■= [0, 1) \ D = |J 



1 



1 



n + 1 n 



Denote Z) = in the ECF case, respectively D = Do in the OCF case. In both ECF or 
OCF situations the signs e n = e n (x) and the digits a n = a n (x) are given, for x G £1, by 

eo = 1, e n = e£>(t n _i), a = 0, a n = a fl (t„_i), 

where t n = t n (x) = T^(x). On the D-continued fraction expansion the iterates of the Gauss 
type map Tp act as a shift map by 

rg[[(ai,ei), (a 2 ,e 2 ), •••]] = [[(a„+i, e n+ i), (a n+2 , e n+2 ), . . •]], Vn G N . 

The D-convergents ^ are defined by 



J P-l = 1, P0 = 0, Pn = anPn-i + e„_ip„_ 2 , 

l^-i = 0, g = 1, <7n = a n (ln-i + e n -2q n -2, 
or in equivalent formulation 



Pn-l Pn 



Pn-2 Pn-l 

'o e \ A) ei N 
1 aW ll a 2 



e n _i 

1 a„. 



(2.2) 



(2.3) 



e n _i 

1 a™ 



Vn G N. 



The following elementary fundamental relations are satisfied: 

Pn-iQn ~ PnQn-i = (-l) fc e ei ■ ■ ■ e„_i =: 5 n , 

Pn-l Pn 5n 



Qn-1 Qn Qn-lQn 
Pn ~\~ Pn—l&ntn 



-, Vn € N , 



x = 



Qn ~i~ Qn—l&ntn 

The latter equation is equivalent to 

&ntn = G n Tjj(x^ = 

Upon (2.4) we infer 

Qn% Pn 



Vn G N. 



Qn% Pn 

-q n -ix+p n -i 



Vn G N. 



< 



< 1, Vr G n, Vn G N. 



(2.4) 



(2.5) 



-q n -ix +p n -i 

It is well-known and plain to check for every continued fraction that if x is as in (2.1), then 

Qn-l 



Qn 



[[(an, e„-i), (a n -i, e„_ 2 ), . . . , (a 2 , ei), (cti, *)]], Vn G N, 



(2.6) 



where (a±,*) means that the finite expansion terminates with a\. 
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In GL2(Z) consider the matrices 

j= 1° ; i . \ i : i . /• i 







and denote their images in »SL2(Z/2Z) by [7], [J], [A], [B]. Clearly {[/], [J]} forms a subgroup 
on two elements of SX2(Z/2Z) and {[/], [A], [B]} forms a subgroup on three elements of 
SX2(Z/2Z). Consider the sets 

He := {M G 71 : 1 s: Q < Q', M = I or J (mod 2)}, 
TZo ■= {M £TZ: X M > g, M = I, A, or 73 (mod 2)} . 
For M G GL 2 (Z) denote 

O' O'r - P' 

Am = |-, E M(x)= _ Qx + p , xiQ. (3.1) 

3.1. Successive convergents for ECF(x). 

Lemma 3.1. In the ECF expansion, ^ qk-i ^ 1> Pfc+i ^ Pfe ^ 1> Qk — Pk ^ 
- Pfc-i ^ 1 /or every k ^ 1. 

Proof. Let (x n ) be a sequence defined by x n = a ra x n _i + e n _ix ra _2 with a n an even positive 
integers and e n G {±1}- Suppose that Xk ^ %k -i ^ 1 f° r some ka ^ 1. Then Xfc 0+1 ^ 
2xfc — Xfe _i ^ Xk - This shows inductively that x n ^ x n _i ^ 1 for every n ^ fco- The 
statement follows by taking (x n ,ko) = (q n ,l), (x n ,ko) = (p n ,2), and respectively (x n ,ko) = 
{q n -PnA)- □ 

Furthermore, since Vn-iQn — VnQn-i = il, it follows that g n (x) > (7 n _i(x), for all n ^ 2 
and x£fi. 

Proposition 3.2. For each x & CI the following are equivalent: 

(i) ^,^7 successive convergents in ECF{x). 

(ii) M = J) G 7^ and < \E M {x)\ < 1. 

Proof, (i) => (ii) Suppose M = (jj) = (^Ii'S) f °r some n ^ 1. From Lemma 3.1, 

(l 6 afc 1 ) = J (mod 2) and equality (2.3) we infer that M G 1Ze- The second condition in 
(ii) follows from (2.5). 

(ii) ==> (i) Consider first the case Q = 1. Only the matrices M = (iq') an d = 
(i ^Q' 1 ) ma ^ ar ^ se - Since M = I or J (mod 2) only the former case can occur and Q' 
is necessarily an even positive integer. The corresponding inequality < | | < 1 is 
equivalent to x G ( q^ +1 , <jt) U (^7, qtzj ) or, according to the definition of 01, to a\ = Q', 
showing that j, ^7 are successive convergents of x. 



6 



FLORIN P. BOCA AND JOSEPH VANDEHEY 



When Q > 1, take (I ^ 1): 

e M = 1, Qo = Q' ~ 2£Q, P Q = P'- 2iP if [A] = 2£, 
e M = -1, Qo = 2£Q - Q', P = 2£P - P' if [A] = 21-1, 



M = 



Po P y 
Qo Q i 



In both cases one has < Qo < Q, M = M ( ° e 2 *f ) , and so M = / or J (mod 2). Since 
Q' > Q > Qo, the condition < \Em\ < 1 is equivalent with x lying between q,^_q and q>Zq j 
while < |-EmoI < 1 is equivalent with x lying between q+q° and qZq° . When |j < |j7 the 
former implies the latter because of 

P - P _ {2£ + 1)P -P' P P' + P P' P' - P 
Q-Qo ~ (2£ + l)Q - Q> < Q < Q' + Q < Q 7 < Q> - Q 
P + Po _ P'-(2£-l)P Pp _ P'-2£P 
Q + Qo~ Q' ~{U- l)Q Q ~ Q' - 21Q wiien [AJ " U ' 

and of 

Po _2£P-P' P + P _ (21 + l)P - P' P P' + P P' 
Qo~ ~ 2£Q - Qi < Q + Qo ~ (21 + l)Q - Q' < Q < Q' + Q < ~Q' 
P'-P P-Po P'-(2£-l)P 
< Q^Q^Q^Q- = Q'-(2£-DQ when W = ^-l. 

When j^r < ^, analogous inequalities show that < \Em\ < 1 implies < \Em \ < 1- 
Furthermore, the inequalities ^ Po ^ P follow from \P'Q - PQ'\ = \PQo — PoQ\ = 1 and 
P ^ 1. □ 

3.2. Successive convergents for OCF(x). Denominators of successive convergents for 
OCF(x) satisfy ([11, Eq. 2.10]) 

Qn 

r n ■= 

Qn-l 

= a n + e n _i[[(o„_i,e n _ 2 ), (a n _ 2 ,e„_ 3 ), . . . , (a 2 ,ei), (ai, *)]] 
> a n - [[(3, -1), (3, -1), . . . , (3, -1), (3, *)]] ( 3 - 2 ) 
>a n - [[(3,-1), (3,-1), (3,-1)...]] 
= a n - 1 + 1/G = a n - 2 + G. 
In the opposite direction one has 

r n = a n + <a n ^ &n ~l ^ a n + — !— - = a n + G. (3.3) 

r n _i a n _i - 2 + G G — 1 

In particular (3.2) and (3.3) show that if a n ^ 3, then r n > 1 + G, proving 

Lemma 3.3. If r n ^ 2 + g then a n = 1, and in particular e n = 1 one? < - q ^x+p n ^i < ^ 

Proposition 3.4. For eac/i x £ Q the following are equivalent: 

(i) |j,|j7 successive convergents in OCF(x). 

(ii) M = |gg;j G T^-o and one of the following two conditions holds: 

(*) X M ■= % > 2 + g and < \E M (x)\ < 1. 
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(**) g < Am ^ 2 + g and < E M (x) < 1. 
Proof, (i) => (ii) Suppose that there is n ^ 1 such that 

M=(° e ° = 1 )(° ei )...(° e n-l\( P n-l Pn\ ^ 

yi ai y yi a 2 y a n y y<?n-i g n y 

Since (? ^T 1 ) = (? i) = A (mod 2) and {[/], [A], [B]} forms a subgroup of SL 2 (Z/2Z), it 
follows that M = I,A, or B (mod 2). The inequality GQ' > Q follows from (3.2), while 
^ P = p n -i ^ Q = q n ~i, < P' = p n ^ Q' = q n are well-known (they follow as a result 
of the RCF — > OCF algorithm or can be directly deduced from p n _iq n — p n q n _\ = ±1). 
Properties (*) and (**) follow from (2.4), (2.5), and from Lemma 3.3. 
(ii) => (i) Consider the partition 

(p, oo) = S\ U S2 U £3, where 

Si = (p,l)U(2 + p,3)U(4 + p,5)U..., 

5a = [1,2) U [3,4) U [5, 6) U..., 

S 3 = [2,2 + p)U[4,4 + p)U[6,6 + p)U... 

For each matrix M = |g^j G TZq with A = Am, define 



kM 



Note that 



21-1 if A€«S 2 , [A] = 2£-l, O 1, 

2^+1 if A G Si, [A] = 2t, O 0, and {A} > g, 

21-1 if A G S 3 , [A] = 2£, £ ^ 1, and {A} < p. 



fc M ^ 3 ^> A > 2 + p = G 2 . 



We prove the following statement: 

Lemma 3.5. Let x G $7 and M = ^gg,j G 7£o Q = mm{Q,Q'} > 1 and satisfying 
(*) or (**). T/iere exisi eM G {±1} and Mo = G ^-o s^c^ 

M = M (; «-), (3,) 

eM+^Af,, ^ 2, Mo satisfies the corresponding property (*) or (**), and Qo = rnin-fQo, Q} ^ Q- 
Furthermore, if A = Am G Si US2, f/ien we can ta&e Qo < Q- 

Proof of Lemma 3.5. Consider the following integers: 



eM 



1 ifAGS 2 US 3 , 
-1 if AG Si. 
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Qo 



Q'-(2£-l)Q if A G <S 3 , [A] 
Q'-(2£-l)Q if A € <S 2 , [A] 
(2£ + l)Q-Q' if A G Si, [A] 

(1 + {A})Q ifAGS 3 , 
{\}Q if A G S 2 , 

(1-{A})Q if A G Si. 



2€, £ > 1, and {A} < g, 

2£-l,£^l, 

2£, £^0, and {A} > g. 



2£, £^1, and {A} < g, 

2£-l,£^l, 

2£, £^0, and {A} > g. 



P'-(2£-l)P if A G S3, [A] 
P'-{2£-\)P if A G 52, [A] 
(2£ + l)P-P' if A G <Si , [A] 
Equality (3.5) holds in all cases with this choice for Qo and Po- One plainly checks that 

(2 + 5,00) if A € Si, 
G < (l,oo) if A G S 2 , 
(5,1] if AG S3. 



Q_ 

Qo 



In particular this shows that Ao > g- The inequality eM + &M ^ 2 is trivial when A G S2US3. 
When A G Si we have Aq > 2 + g hence Jcmq ^ 3 and eM + &m ^ 2. 



Clearly 



e M 

1 &M 



= A (mod 2). The inequalities ^ Po ^ Qo follow immediately from 

PoQ — PQo = ±1 and P < Q, the latter one being a consequence of the assumption Q > 1. 
The fact that Mo satisfies either (*) or (**) follows from Lemma 3.6. □ 

Back to the proof of Proposition 3.4, note that when A G (g, 1] one has < Qo = Q — Q' < 
Q' < Q (the first inequality holds because G < 2), while for A G (Si U S2) \ (g, 1) it is plain 
that < Qo < Q < Q' ■ Hence whenever A G Si U S2 one has min{Qo) Q} < niinjQ, Q'}- 

When A G S3 one only has min{Qo,Q} = minlQ,^'} (actually Q < Qo < Q')- However, 
in this case eM = —1 so kM ^ 3, and Xm = <j- G (g, 1). Thus one can apply the same 



procedure to Mq and find M_i 



P-i Po 
Q-i Qo 



G IZq that satisfies (*) or (**), and such that 



eM + ^m_i ^ 2, and Q-i := min{Q_i,Q } < Qo = Q (this inequality 
is strict because Ao G (g, 1) C Si). 

We next discuss the case Q = 1. When Q' = 1 ^ Q, the inequality = ^ ^ 9 yields 
Q = 1. Hence M = (§ 1), with y, y successive convergents of every x G (0, 1) that satisfies 
< ^ < 1, i.e. of every x G (\, l). Suppose now Q = 1 <Q'. When £j<^ = Q'<2 + 0, 
one has Q' = 2 and only the matrices M = (12) an d M = ( { \ ) may arise. But the 
former matrix is not admissible being = ( \ q) (mod 2), while the latter matrix corresponds 
to < < 1, hence x G (^, |), ei = 1 and a\ = [^] = 1, and indeed \, \ are successive 

convergents in OCF(x) for every x G (5, |). When 2 + g < ^ = Q' the only matrices that 
may arise are M = ( ? q' ) with Q' ^ 3 odd, and respectively M = ^ j ^ with Q' ^ 4 even. 
The inequality for the former is < | < 1, which gives x G ( q} +1 , 757) U (757, qtzi ) 



with Q' odd, so that a± = Q' (and ei = 1 respectively e± = —1). The inequality for the 
latter is < | ^'l'^'* 1 1 < 1, giving > x > ^rf ^ | so ei = 1, a\ = 1. Furthermore 

one has ^ < ^ — 1 = Td(x) < qr^ with Q' — 1 ^ 3 odd integer, so 02 = Q' — 1 and 
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M = (1 1) (1 q'-i)i showing that indeed y, Qq, 1 are successive convergents in OCF(x) for 
every x with q7 +1 > x > q, _ 1 and Q ^ 4 even. 

This inductive process on Q now implies that (3.4) holds for some e±, . . . , e ra _i G {±1} 
and ai,...,a n odd positive integers with a + ^ 2, Vz G {1, . . . , n — 1}. Conditions (*) and 
(**) show that x lies between ^"~^ n ~ 1 and ^"t!?"' 1 when q n > o n _i, and between ^ and 

^iXqlZl wnen ?n < So x is of the form [[(ai,ei), (a 2 ,e 2 ), . . . , (a n _i,e n _i), (a n + i, *)]] 

for some t G (—1, 1) when q n > q n -i, and respectively t G (0, 1) when q n < q n -\. Therefore 
tt^ 1 = r>, 7T = 7V are successive convergents of x. □ 

<2n-l Q' fc v ° 

Lemma 3.6. W^i/i the definitions from the proof of implication (ii) =>■ (i) in Proposition 
3.4 one /ias 

(i) If g < A < 1 iften < < 1 \E Mo (x)\ < 1. 

(ii) // 1 < A < 2 + 5, then < £ M (x) < 1 < £ Mo (z) < 1- 

(iii) If2£ + g < A < 2£ + 1, £ ^ 1, iften |^m(^)| < 1 => -1< E Mq (x) < 0. 

(iv) //2£ - 1 < A < 21 + g, i > 2, then \E M {x)\ < 1 => < E Mo (x) < 1. 

Proof. In all cases < Em{x) = ®qx+p < 1 is equivalent with x lying between ^7 and g;+g ; 
while < Em (x) = Jq^+p^ < 1 is equivalent to x lying between g and q+q° ■ 

(i) In this case Qo = Q — Q' < Q and — 1 < Em (x) < 1 is equivalent to x lying between 
= MfS and |fq\ = W ■ The cond usion follows because fgfg < g < g^g < g 

when € < 5 and § < < € < 1f5 when § <§■ 

(ii) In this case q+q" = Q" 7 ana - x between gy and gr{-g implies x between g and gy. 

(iii) In this case < Qo = (21 + 1)Q - Q' < Q < Q' , — 1 < Em(x) < 1 is equivalent to x 
lying between g/+g and q/Zq j and — 1 < Em (x) < is equivalent to x lying between g and 

f=?£ = W^m- The implication follows because either g < g^g < g < g^g < g^Ig 

P'-2£P , P'-P EL ^ P'+P E. 
or Q'-2*Q Q'-Q ^ Q' ^ Q'+Q ^ Q" 

(iv) In this case Q' > Q and q+q° = Q'Z^t-2)Q ■ ^ ne implication follows because — 1 < 
Em{x) < 1 is equivalent with x lying between g,^g and g,~g , < Em {x) < 1 is equivalent 

1 P 1 P-t-Pn l-l P P'-t-P P' P'—P P' — (2£ — 2)P 

with x lying between ^ and and either ^ < ^ < w < q^q < q/-(^- 2 )q or 

P'-m-2)P P'-P p^ P'+P P , p, 
Q'-(2£-2)Q ^ Q'-Q ^ Q> ^ Q'+Q ^ Q n °lQb. LJ 

The following statement will also be useful: 

Lemma 3.7. Denominators of successive convergents in OCF satisfy 

(i) Qn+2 > q n - 

(ii) q n+3 > q n . 

(iii) q n+ 2 > mm{q n ,q n+ i}. 

Proof. By Proposition 3.4 and its proof ^n+i > 2 > 2o > 1, G (1, 2) =^> > 

1 ^> ^±2 > 1 and £n±2 G / ft ^ £n±i > 2 + n ^ «n±2 > 0(2 + o) > 1. Thus in all 
possible cases g ra+ 2 > q n , which establishes (i). 

(ii) follows from ^g(#,1)^^>2+ 5 ^^>(2+ g)g 2 = 1, 2=±s = A € 

V ' ; <?n+l A-l Q„ A-l» y ' <J„+2 V ' ilJ qn + 1 Vi " ^ 

«n+i>2 + o^^±3> 2o(2 + o) > 1 and 2e±s > 2 + o 2a±3 > (2 + g)g 2 = 1. 
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To prove (iii) suppose that g n + 2 ^ Qn+i- Then ^±2. g x), which gives in turn > 
2 + g, and therefore 2=±a > g{2 + g) > 1. □ 

Remark. Proposition 3.2 was originally proved, using a different method, by Kraaikamp 
and Lopes [7], but Proposition 3.4 is, to the best of our research, new. Our proofs have an 
additional benefit of implying how to derive a n and e„_i (and hence q n -2) if only q n -i and 
q n are known. 

Our investigations yielded yet another method of proof, significantly longer but more direct, 
which we sketch here. Examples 1.8 in [8] explain how to algorithmically generate the OCF 
expansion of x from the RCF expansion of x using insertion 

(replacing [[..., (a n , 1), {a n+1 , e n+1 ), . . .]} 

with [[..., (a n + 1, -1), (1, 1), (a n+ i - 1, e n+2 ), . . .]]) 

and singularization 

(replacing [[..., (a n , e n ), (1, 1), (a n+2 , e n+2 ), . . .]] 

with [[..., (a n + e„, -e„), (a n+2 + 1, e n+2 ), . . .]]). 

Both of these operations alter the sequence of convergents: insertion adds a new convergent, 
while singularization deletes one. Nevertheless, it can be shown that if ^, ^ are successive 
RCF convergents to some x, then either are successive OCF convergents to x or 

® q/ P ar e successive OCF convergents to 1 — x. (Only one of these pairs forms a 
matrix that is congruent to /, A, or B modulo 2.) By carefully following how insertion and 
singularization change the last e n _i and a n in the RCF expansion of ^7 into the last 6m— 1 
and a m of the OCF expansion of £7, we can determine exactly what e(M) and a{M) must 
be and hence how to derive Pq and Qo- A similar proof works for the ECF case as well. 



4. Estimating the limiting joint distribution for ECF and OCF 

For each M = ^ ^ ) e K and f € (0, 1] denote by l£{M) (respectively, I^(M)) the set 
of solutions x of ^ Em(x) ^ £ (respectively, of — £ ^ Em{x) ^ 0). The Lebesgue measure 
of lf{M) is 



ft(Q,Q') 



P' ± £P P' 



The integral 



Q'(Q'±£Q)' 

min{x3i;,a;iR} 



/■oo fmm(3;3tj,a;iJt) 

F± = F±(xi,x 2 ,x 3 ,X4) := dv duf±(u,v) 

JR/xo Jo 



IR/X2 JO 



v ± X4 min{x3V, Xii?} 



ttJ ± X3X4 min{w;, x\ } 



can be expressed when X3 ^ xix 2 as 

fXiX2XA ^ 



J 

Jo 



log(l ±t) = T Li 2 (T^i^2X4), 
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and when X3 < X1X2 as 



f dw f dw w± £1x3x4 
F±= — log 1 ± X3X4 ± / — log 

J X3 /X2 W J X1 W W 

= ± log(l ± X3X4) log =F Li 2 (=Fx 3 x 4 ), 

so F± is as in (1.5). 

4.1. The ECF case. By Lemma 3.1 and Proposition 3.2, for each R > 1 and x G f2 there is 
a unique M = ^g^j € 7£g with Q ^ R < Q' and |£m(x)| < 1. Given xi, X2, X3, X4 G (0, 1) 

consider A^'^^^lx, R), the number of matrices M G 7£e that satisfy (1.1) and (1.2). One 
has 

£ ' (F) = £-xi,X2,X3,X4,{R) = / ^X\,X2,Xz,Xi i X l ^) dx . 

Jo 

For T G {I, J, A, B} we shall estimate 

£±(R):= £ fxMQ')- 

Q'^R/x 2 
Q^min{x3Q' ,xiR} 
M=T (mod 2) 

This can be done by Mobius summation, as in the following standard lemmas (for Lemma 
4.2 see, e.g., [2, Lemma 2.1]). 

Lemma 4.1. For every interval J, every function g G C l (J) of total variation Tjg, and 
every integer x, with ctq the divisor counting function, 

£ 9{a)= £ g(a) = ^ f g(u)du + O(a (q)(\\9\\oo + Tjg)). 
aeJ,be[i,q] aeJ q Jj 

ab=x (mod q) (a,q)=l 
(a,q)=l 

Lemma 4.2. For every interval J, every V G C 1 [0, N], and every £ G N 

£ — %) = cw rn^^ + o^dlFiu + To^iogiv), 

K^jv q Jo 
(q,e)=i 

with 

-1 



per 

p\t 

Changing b to q — b in Lemma 4.1, we infer 

Corollary 4.3. Suppose q is an odd positive integer. For every interval J, every g G C 1 (J) ; 
and every integer x 

£ 9(a) = — I g(u)du + O((\\g\\oo+Tjg)a (q)). 

aeJ,be[l,q/2] Q 3 

ab=x or —x (mod q) 
(a,q)=l 
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Since P'Q — PQ' = ±1, P',Q even and Q' odd entail P odd, we infer (with Q = 2q, 
P' = 2p', x the multiplicative inverse of x (mod Q')) 



£f( R )= E E ftS^Q') 

Q'^R/x 2 qe[l,min{x 3 Q',yiR}/2] 
Q'=l (mod 2) p'€[L,Q'/2] 

p'q=±A (mod Q') 

y\ /•min{x 3 Q',xiR}/2 



E l^iri f^Q')d q +o £ (Q'-^)\ 

(4.1) 

1 f(Q') f min i x sQ' \ x i R } 



Q'^R/x 2 
Q'=l (mod 2) 



E ^ 4 (.,Q0^ + O £ (^ 1+£ ) 



2 

Q'=l (mod 2) 



On the other hand, P'Q—PQ' = ±1 and Q' even entail that both Q and P' are odd, and the 
condition P even is equivalent to P'Q = ±1 (mod 2Q'). Since in this case (p{2Q') = 2tp(Q'), 
we infer 



£j(*)= E E fx 4 (Q,Q') 

Q'^R/x 2 Qe[l,mm{x;iQ' , Xl R}] 
Q'=0(mod2) P'e[l,Q'] 

P'Q=±1 (mod 2Q') 



^ / <jj(20') r^H x 3Q',xiR} 

E ( ^ y o /£(«■ tf) d - + o £ (<r 2+£ ) 



Q'^R/x 2 
Q'=0 (mod 2) 



leading to 



C E >±(R) = Cf(R) +£±(R) = ^± + £ (R-^), 



and concluding the proof of (1.3). 

The corresponding estimates for C^(R) and C^(R) are useful for the OCF situation. To 
estimate Cg(R), note that P'Q — PQ' = ±1 and Q' even entail that both P' and Q are odd, 
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<p{2Q') = 2ip(Q'), and thus 

E E fitSQiQ') 

Q'^R/x 2 Qe[l,mm{x3Q',xiR}] 
Q'=0 (mod 2) P'e[l,Q'],P'Q=±l (mod Q') 

^1=1 (mod 2) 

El E ftMQ')- E 



Q'^R/x 2 \ Qe[l,min{x 3 Q',xiil}] Qe[l,min{:r3Q',:riP}] 
Q'=0 (mod 2) P'e[l,Q'],P'Q=±l (mod Q') P'e[l,Q'], P'Q=±1 (mod 2Q') 

r<min{a;3Q' xiP} 

E KT-ir)/. /£(»,Q')A l + o,(Q'- 2 -. 

»'=0 (mod 2) 



E 



Q'^R/x 2 
Q'=0 (mod 2) 




(4.2) 



(cM- c,(2) ) F± + ° e( ^ 1+£) 



3C(2) 

Finally, P'Q — PQ' = ±1 and P even entail that both P 1 and Q are odd, and so 

Ca( r )= E E fxMQ') 

Q'^R/x 2 Qe[l,min{x3<3',xiJi}] 
Q'=l (mod 2) P'e[l,Q'],P'Q=±l (mod 2Q') 

( mOO'\ f m in{x 3 Q' ,X!R} \ 

E #/ fi(u,Q')du + O e (Q'-^)) (4.3) 

Q'=l (mod 2) 

= ^ F ± + O e (/2- 1+e ) = + O e (^ 1+£ )- 

4.2. The 0(7F case. This requires more caution as the sequence of denominators of suc- 
cessive convergents is not monotonically increasing in general. We wish to characterize those 
matrices M G IZo for which ^, ^7 are successive convergents of x G Q and Q = 9n H ^ R < 
Q' = (MjH-i- A priori, Lemma 3.7 shows that for each R > 1 there is at least one pair and 
at most two pairs (Q,Q r ) of denominators of successive convergents of x with Q ^ R < Q' . 
Moreover, if there are two such pairs (Q,Q r ), then they must be of the form (q nR ,q nR+ i) or 
{qn R +2, Qn R +3)- We wish to precisely distinguish ur from n# + 2. Because all predecessors 
of Qo in the sequence of denominators of OCF convergents are < Q by Lemma 3.7, equality 
(Q,Q') = (q nR ,q nR+ i) occurs exactly when 

Q ^ R<Q' and R > Q . 

Note that if A = ^- G S1US2, then necessarily Q > Qq. Furthermore, if A G S3, then Q < Qq. 
The contribution of those pairs (Q,Q') with A G £3 and Qq = Q(l + {A}) > R should be 
subtracted, and so we can write 

£°>+{R) = C+{R) + C+(R) + C+{R) - V^R), 

£°>-{R) = CJ(R) + C A (R) + C B (R) - V 2 (R) - V 3 (R), 
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with 

X4 



t>i{R)= E fl(Q,Q') = Y, E 



Q'(Q' + x 4 Q) ' 

M£K ,Q'>R/x 2 l>\ M€Ro,Q'>R/x 2 y ' 

Q^min{x3Q' ,xiR} Q^min{a;3Q',a;i_R} 
X=Q'/QeS 3 , Q(1+{X})>R 2£Q^Q'<(2£+g)Q 

Q'>R+(2£-l)Q 



Q^min{x3Q' ,xiR} 
A=Q7Qe[2,2+g), Q'>R+Q 



v-s(R)= E = E E 



X4 



Men ,Q'>R/x 2 02 M&n ,Q'>R/x 2 K ' 

Q^min{a;3Q',a;i_R} Q^min{a;3Q',a;ii?} 
A=Q'/Qe<S 3 ,A>G 2 2£Q^Q'<(2£+g)Q 
Q{1+{\})>R Q'>R+(2£-l)Q 

Clearly T>2(R) = when min{:ri:r2, £3} ^ g 2 . When min{2;iX2, X3} > g 2 , the method 
employed in (4.1), (4.2) and (4.3) applies, leading to 

V 2 (R) = £^^^ + £ (R-^), 

with D2 as in (1.6). 

The estimation of T>i(R) is slightly more involved because £ can take infinitely many values. 
Note that T>i(R) = unless min{xiX2, X3} > 2T+g- eacn £ N consider the integral 



inm =- J/ 



X4 du dv 



v(v + x^u) 

v^R/x2,u^mhi{x:jV,xiR} 
2£u^v^(2£+g)u 
v>R+(2£~l)u 

The change of variables (v,u) = (Ry,Rx) shows that I^(R) does not depend on R and is 
given by (1.7). Note also that 

/•xi A2£+l)x , , 

lj(R) < / dx / -| « -. (4.4) 

Jo J2lx V *■ 

A trivial estimate yields 

E E E E 



Q'(Q' + x A Q) ^ Q 2 £ 2 

R/x 2 ^Q%(2£+l)R 21+1 21 1^Q%(2£+1)R 21+1 21 

El \ - x - 1 logi? 

72 E E 02« 



£2 Q2 ^ pl/2 ' 

l^R 1 / 2 Q£{l,2R]Q>£[2£Q,{2£+l)Q]^ 
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and thus in the definition of T>±(R) we may take I € [1,-R 1 / 2 ] inserting an error term <C 
RT 1 ! 2 logi?. Employing Lemma 4.1, the resulting main term can be expressed as 

^ ^ U' k Q'(Q' + x A u) + ^ } j 

KB}/ 2 Q'^R/x 2 \ ^ ^ 4 ' J 



Q'<(2£+g)x 1 R 



>£<i?V2 Q'^R/x 2 21+g 
Q'<(2e+g)x 1 R 



Employing now Lemma 4.2, the main term above becomes 



and so 



V ^ = ck) E ^ + + £ ( J R- 1/2+£ ). (4.5) 
From (4.5) and (4.4) we eventually infer 

The sum T>s(R) is similarly estimated as in formulas (1.6) and (1.7). 

5. Joint distribution for Nakada's q-expansions 

We illustrate how explicit renewal type results can be obtained in the case of Nakada's 
a-expansions NCF a , a € \_\, l] . Such continued fractions, defined in [10], have been studied 
in [10, 6]. Here the unit interval is replaced by Q a = [a — l,a) and the Gauss shift by the 
map T a : Q a — > f2 a defined for x ^ by 3 

T a (x) 

A construction of the natural extension T a on a space Q, a C M 2 , together with an explicit 
invariant Borel probability measure fi a on [l a was found by Nakada [10]. He also proved that 
(QctiTa, Ha) is a Kolmogorov automorphism. With g = ^ = 1 — g 2 the set is given for 
g < a 1 by 

[a- 1,(1 -a)/a\ x [0,1/2) U ((l-a)/a,a) x [0,1] U [a- 1,0) x {1/2}, 
and for | ^ a ^ g by 

[a - 1, (1 - 2a) /a] x [0, g 2 ) U ((1 - 2a) /a, (2a - 1)/(1 - a)] x [0,1/2) 
U ((2a - 1)/(1 - a), a) x [0, g) U [- 5 2 , (1 - 2a) /a] x {<? 2 } 
U ((1 - 2a)/a,0) x {1/2}. 



1 






1 


+ 1 - a 










X 











'Here we use the notation from Sections 5 and 6 of [6] . 
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Kraaikamp's thoughtful analysis (see especially Theorem (5.3) and Definitions (5.7) and (5.8) 
of [6]) also provides characterizations of pairs of successive convergents for such continued 
fractions if a G [|, l] . 

Proposition 5.1. For each x G Q a \ Q the following are equivalent: 

(i) ^, ^ successive convergents in NCF a (x) with Q,Q' > 0. 

(ii) M=(2£) GGL 2 (Z) and (^ M (x),l/A M ) GO*- 

This dynamical system was studied by Kraaikamp [6] in the more general setting of S- 
expansions, and the above proposition can be likewise generalized if we replace NCF a (x) with 
CFs(x), the S-expansion of x, and replace £l a with Q s , the space of the natural extension 
associated to S. 

We wish to estimate the Lebesgue measure >C Xl j X2)X3lX4 (-R) of the set of numbers x G Q a \Q 
for which there exist ^, ^7 successive convergents in NCF a (x) that satisfy (1.1) and (1.2). 
We shall require that xi,X2,x$ are in the set (0, 1] if g < a ^ 1, n the set (0, 1/2] if a = g, 
and in the set (0,g] if 1/2 < a < g; moreover, we require X4 G (0, a] when we look at £ + 
and X4 G (0, 1 — a] when we look at Cr . The set is a union of rectangles and horizontal 
line segments, but we may ignore the line segments for large R: in particular, the inequality 
Q' ^ j| shows that the pair (Q 1 , Q) = (2, 1) makes no contribution to £ ± for R > 2, so the 
situation A M X = \ can be ignored, and Am is always rational, so the situation A M X = g 2 can 
also be ignored. As a result, the cases that appear in £x°^ )X3jX4 (.R) for i? > 2 are exactly 
the following: 



For 5 < a ^ 1: 



For 1/2 g: < 



Am = Q'/Q > 2 and a — 1 ^ Em(x) < a, or 

1 < Am < 2 and ±=2 < E M (x) < a. 

Am > G 2 and a — 1 ^ Em{x) < a, or 

2 < Am < G 2 and < £m(x) < a, or 
G < Am < 2 and < E M {x) < a. 

The varying lower bounds on Am depending on the value of a are the reason for our case-based 
restrictions on the values of xi,X2,x^. 

Let £+ X2 X3 X4 (a; R) denote the Lebesgue measure of the set of numbers x G [0, 1] \ Q 
for which there exists M = (q q>) G GL 2 (Z) with Q,Q' > 0, ^, ^ G [a - l,a) and (1.1) 
together with ^ ~q x +p ^ x 4 hold. The corresponding set where the latter inequality is 
replaced by —X4 ^ ®q~+ p ^ is denoted by £~ lX2X3X4 (a;R). In both cases, x\, X2, x%, X4 
are parameters in (0,1]. When a = 1, it is clear that C + is exactly the joint distribution 
considered in [16] (where the notation used is N(R)). However, by the following equation 

± < . m _ \ ^ X4 



£-x 1 ,x 2 ,x 3 ,x 4 ( a i R) Z^Z 



Q'(Q' ± X4Q) 

Q'^R/x 2 Qe(0,mm{x 3 Q',x 1 R}} ^ v ^ ^' 
P'e(a-1)Q'+[0,Q') 
P'Q=±l (mod Q') 



f)t(f)> 4- T ./11 ^Xl,X2,X3,X4(-^)' 



Q'(Q' ± x 4 Q) 

Q'^R/x 2 Qe(0,min{x 3 Q',xi.R}] ^ y ^ ^ ' 

(Q,Q')=i 
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we see that £ ± (a; R) does not depend on a. As R tends to infinity, converges to 2F ± /C(2). 

The joint distributions £( a )' ± and can now be directly related as below. For the sake 
of space and readability we omit the appearance of x\, X2, and R, which are assumed to be 
the same on the left- and right-hand sides of the equations. 

When g < a ^ 1, we have 

r(a),+ = J ^niin{a;3,l/2},X4 if ^ £4 ^ 

[^xs^ rx 3 ,(l-a)/a + ^min{x 3 ,l/2},(l-a)/a 11 a ^ ^4 < a, 
4*5x4 = £ min{x 3 ,l/2},X4 if < X ± < 1 ~ 

When i ^ a ^ g, we have 



^min{x 3 ,l/2>,X4 if ^ X 4 ^ 2 ™_^ , 



/•(<*)>+ _ 

^£3, £4 A 



^ (<*)>- 



"'min{3;3,l/2},2:4 

r+ _ r+ 

*-x 3 ,X4 t, x 3 ,(2a-l)/(l-o;) 

+^min{x 3 ,l/2},(2a-l)/(l-a) if l^T ^ ^4 < «• 

£~ , im if < x 4 < 2s=I 

mm{x3,l/2},X4 \ <± \ a 1 

r- 4. £- 

min{x3,3 2 },X4 ^ min{x3,l/2},(2a— l)/a 

— £~. r wo in / if 2(1-1 ^ Xa < 1 — a. 

mm{x3,<?' ! },(2a— l)/a a ^ * ^ 



Recall that X3 ^ g in this case. 
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